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In the usual Achlioptas processes the smallest clusters of a few randomly chosen ones are selected 
to merge together at each step. The resulting aggregation process leads to the delayed birth of 
a giant cluster and the so-called explosive percolation transition showing a set of anomalous fea¬ 
tures. We explore a process with the opposite selection rule, in which the biggest clusters of the 
randomly chosen ones merge together. We develop a theory of this kind of percolation based on the 
Smoluchowski equation, find the percolation threshold, and describe the scaling properties of this 
continuous transition, namely, the critical exponents and amplitudes, and scaling functions. We 
show that, qualitatively, this transition is similar to the ordinary percolation one, though occurring 
in less connected systems. 
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Aggregation processes based on progressive merging 
together the smallest clusters of a few randomly chosen 
(Achlioptas rule) have attracted much attention in the 
last years M- These specific processes and models 
show a number of unusual features E3 distinguishing 
them sharply from ordinary aggregation processes and 
standard percolation, in which random clusters merge 
together with probability proportional to their sizes [18l — 
;20]. Apart from the delayed percolation phase transition, 
which is continuous, as we found for a wide range of sys¬ 
tems [d, ; 9ij and which was proven mathematically fl2l . [l3j . 
these models demonstrate a uniquely small exponent /3 
of the percolation cluster and unusual scaling functions. 
The smallness of /? makes the transition so “sharp” that 
it is difficult to distinguish it from discontinuous in sim¬ 
ulations, which resulted in the term “explosive percola¬ 
tion” [l|. A few real-world applications of these processes 
were identified !2ll - l23l |. The Achlioptas processes, gener¬ 
alizing percolation, constitute a wide class including the 
processes generated by the original “product rule” (two 
clusters with the smallest product of sizes are selected) 
[ll], the sum rule (clusters with the smallest sum of sizes 
are selected), the rule selecting the smallest clusters 0], 
and many others, of which only a small number were ex¬ 
plored. The problem is how far from the standard perco¬ 
lation scenario can these diverse rules and their variations 
lead? How easy can one deviate from the typical perco¬ 
lation behavior by exploiting the “power of choice” [24} 
in these processes? Notably, in these rules another kind 
of optimization can be considered, namely, selecting not 
the smallest but the largest clusters. In particular, the 
question is: what will happen if we invert the Achlioptas 
rule, that is, at each step merge together the two largest 
clusters of a few randomly selected ones [H, HU? 

In the present article we answer to this question by 
considering a representative set of processes based on the 
inverse Achlioptas rule, for which we derive the Smolu¬ 
chowski equation. By solving these equations numeri¬ 
cally and analytically we find that this rule results in a 
percolation transition taking place at an earlier stage of 


the process, but with the same set of critical exponents 
as in ordinary percolation. We calculate the critical am¬ 
plitudes for the relative size S of the percolation cluster 
and for the size distribution of finite clusters and obtain 
the scaling functions. 

The paper is organized as follows. In Sec. U we intro¬ 
duce the model and describe the evolution equations. In 
Sec. HJ we obtain the critical singularity of the percola¬ 
tion cluster S by using the generating function technique. 
In Sec. m we find the scaling functions and the critical 
exponent r, P(s,t c ) ~ s 1-T . Section [TV] describes the 
order parameter and generalized susceptibility for these 
phase transitions. Finally in Sec.|V]we obtain analytical 
estimates for the transition point and critical amplitude. 
Table U demonstrates a good agreement between the re¬ 
sults of the numerical solution of evolution equations and 
these estimates. 


I. THE MODEL 

We consider the following model incorporating the in¬ 
verse Achlioptas rule and convenient for treatment. We 
start from N isolated nodes. At each time step a new link 
connecting two nodes is added to the network as follows. 
At each step sample two times: (i) choose m > 1 nodes 
uniformly at random and compare the clusters to which 
these nodes belong; select the node within the largest of 
these clusters; (ii) similarly choose the second sample of 
m nodes and, again, as in (i), select the node belonging to 
the largest of the m clusters; (iii) add a link between the 
two selected nodes thus merging the two largest clusters. 
Note that the only difference from our previous works 
is that instead of selecting the two smallest clusters for 
merging BEE na, here we select the two largest. The 
probability distribution P(s,t), i.e., the probability that 
a uniformly randomly chosen node belongs to a cluster of 
size s at time t, gives the complete description of the evo¬ 
lution of this system. Time t is the ratio of the number 
of steps (links) and the number of nodes. We emphasize 
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that, instead of the product rule [lj, we select for merging 
the largest cluster from each of the two sets of m clusters, 
which makes our problem treatable analytically. 

For infinite N, this aggregation process is described by 
the following evolution equation: 


dP(s, 4) 
dt 


s -1 

= s E Q( u > t)Q{s — u,t) - 2 sQ(s, 4), 

14=1 


( 1 ) 


which is the Smoluchowski equation for this process [28j, 
[29!]. Here Q(s, 4) is the probability that a cluster selected 
to merge is of size s. The distribution Q{s, 4) is expressed 
in terms of P(s, 4) as 


Q(s)= 




= mP(s) 


E p ( u ) 

, 14=1 
m —1 


E p ( y ) 


(2) 


Notice the normalization conditions for these distribu¬ 
tions, E s P( s ) = 1 — S and E s Q( s ) = (1— S) m , where S 
is the relative size of the percolation cluster. For large s 
we have 


Q(s) - mP{s) (1 - S)™- 1 , (3) 

both above and below 4 C . Equation m together with 
relation © describe the process exactly in the full range 
of 4, from 0 to infinity. 


II. PERCOLATION CLUSTER SIZE 

Let us define the generating functions 

P(z,t) = EX p (M) (4) 

S 

and 

o-(2,f) = ^2,z 8 Q{s,t). (5) 

S 

We multiply both sides of Eq. CD by z s and sum over s, 
which gives 

9p(z,t) or _,„ ^ „da(z,t) 

~dT = 2 MZ,t)-l)-g (6) 

Using Eq. ©, we find the relation between the functions 

a(z) and p(z) for 2 close to 1, 

cO) = E Q( s ) + E Q( s )( zS ~ x ) 

S S 

s (1 - S) m + m( 1 - S')" 1 - 1 E p (s)(z s - 1) 

S 

= m(l-Sr- 1 p( Z )~(m-l)(l-S) m , (7) 


and so 
dp 


2 [m(l - S) m ~ l P (m - 1)(1 - - 1] 


xm(l — 5)’ 


»-i 


ainz' 


( 8 ) 


By applying hodograph transformation [29j to this par¬ 
tial differential equation we get the ordinary differential 
equation 


din 2 


dt 


SS 2m(l - S) r 


x [(m - 1)(1 - S) m + 1 - m( 1 - S')" 1 ” 1 /?] ,(9) 

which leads to 

In 2 = 2m(m—1) f dt (1-S) 2m_1 + 2m f d4 (l-S 1 )™ -1 

Jt c 

2m 2 pfdt{l-S) 2m - 2 +g(p ), (10) 

Jt c 

where the function g(p) is the initial condition for 
Eq. ©, which can be found from the critical distribution 
P(s,t c ) = aoS 1_T . At 4 = t c , proceeding as in 
for 2 close to 1, we find the singularity 


1 - p(z,t c ) = -a 0 r(2 - f) (-ln;z) T 2 . (11) 

Inverting the function p(z, 4 C ) in this equation we obtain 
_(l_ p )lV(r-2) 


9{p) = - 


_a 0 r(2-r)_ 


which we substitute into Eq. (Hhd- Finally, we set 2=1 
in the resulting equation and obtain 

0 = 2to(to—1) f dt (1 - S') 2 ™" 1 + 2m [ d4 (1 - S)™- 1 
Jtr 


—2to 2 (1 — S) dt (1—S’) 2 


-S 


a 0 r(2-f) 


U/U-2) 

■ ( 12 ) 


From this equation we find the critical singularity of S , 
S “ [-o 0 r(2 - r)] 1/(3 “ r) [2?n 2 (4 - 4 c )] (r - 2)/(3 " r) . (13) 


III. SCALING PROPERTIES 

In this section we find the scaling form of P(s 1 4) near 
4 C using the approach of our previous works [13, 1231 • The 
form of the scaling function is determined by the critical 
exponent r. The scaling function must decay faster than 
any power law, and must take only positive values. We 
show that these conditions are satisfied only for r = 5/2, 
which enables us to find the scaling function for each m. 

Let us obtain the Taylor expansion of P(s,t), 

P(s,t) = Aq(s) + A\ (s)(4 — 4 C ) + A 2 (s)(4 — 4 C ) 2 -I-... (14) 
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(a) 



(b) 



FIG. 1. Functions f(x) of Eq. (1191) for a o = l/\/27r and different values of r. (a) Only r = 5/2 provides a rapidly decaying 
positive function, (b) In the phase below t c (i.e., x < 0), the values r > 5/2 result in functions decaying in a power-law way. 


in terms of r, by sequentially differentiating the evolu¬ 
tion equation m and the relation ([3]) at t = t c with re¬ 
spect to t. Recall that the critical distribution Q(s,t c ) = 
mP(s , tc) for large s, and so the derivatives are 


d n Q(s,t) 

dt n 


d n P(s,t) 


dt 71 


(15) 


Differentiating both sides of Eq. m n — 1 times and re¬ 
placing the right-hand side with Eq. m we find the 
asymptotics of the coefficient A n (s) = d[ n ^P{s,t) |t c /n!. 
Due to Eq. m the equations for A n for each to > 1 
differ only by the factor m 2n on the right-hand side. So 
the asymptotics A„(s) is 


A„(s) s* a n m 2 V- T+ "< 3 - T \ 

with the prefactor a n that we have calculated in 
ordinary percolation ( in = 1), 

2 n [a 0 r(2 - r )] n+1 
an ~ (n + l)!r[(n + l)( 2 -r)]' 


(16) 
13 for 


(17) 


The scaling form of the distribution P(s,t ) is 
P(s, t) = s 1_T ^2 a n [m 2 s 3 ~ T (t - t c )] n 

n 

“s 1 - r /[a 3_T (t-*c)]» (18) 

where f(x) is the series 


f(x) S a 0 T(2 — t ) ^2 

n =0 


[2aom 2 T(2 — t ) x \ " 

(■n + 1)! T[(n + 1)(2 - r)]' 


(19) 


The parameter m appears only as a factor of x. The 
function f{x) depends essentially on r. For x 1 this 
function approaches 0 exponentially, staying positive, for 


any 2 < r < 3. In the phase t < t Cl i.e., x < 0, only one 
value of the exponent r results in a scaling function f(x) 
with the proper decay to 0 as a; approaches — oo. For r < 
5/2, the function f{x), Eq. (fTill) . oscillates around 0 in 
the region x < 0, see Fig. Oja). Since the scaling function 
cannot take negative values, we exclude the range r < 
5/2 from the possible values of r. For r > 5/2, the 
function f(x) stays positive but approaches 0 as a power- 
law as x — >■ —oo, see Fig. [lib). The scaling function must 
decay more rapidly than any power law for x —1 ±oo, and 
so we also exclude the range r > 5/2. 

At t = 5/2 the function f(x) takes the form 

[ 2l n 

—47r ( anm 2 x ) 

- , J 

"—' n! 

n —0 


= ao exp 


—47t (ao m 2 xY 


( 20 ) 


that is, decays exponentially with x 2 on the both sides 
of the transition. Thus r = 5/2 for all to > 1, and 
so Eq. (1201) gives the scaling function of this transition. 
Figure [2](a) shows that near t c at large s the numerical 
solution P(s,t ) of Eqs. ([Tj) and © agrees completely 
with the scaling functions (BUI) . Inserting r = 5/2 into 
Eq. (fl3l) we find 


S = 8nm 2 al(t — t c ) (21) 

near t c . Figure ©b) presents the evolution of the relative 
size of the percolation cluster for each to, which we found 
numerically from Eqs. ([1]) and ©, S(t) = 1 — ^)- 

The curves S(t, m) in Fig.©b) intersect with each other, 
and so Sit) grows slower for larger to above t c . Recall 
that in our model the percolation cluster can be selected 
more than once at the same step. This corresponds to 
adding a new link between two nodes in the percolation 
cluster, which does not changes cluster sizes and happens 
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with probability [1 — (1 — S')" * 1 ] 2 . This probability grows 
rapidly with to effectively delaying the aggregation pro¬ 
cess above t c compared with m = 1. If we forbid the 
same cluster from being selected more than once at each 
step, the delay effect disappears and a larger m results 
in a faster growth of S , which approaches S = t when 
m —> oo. In Table U we show precise results for r, t c , 
and ao, which are computed from P(s < 10 5 ,t) using 
the method of Ref. [§j]. The numerical results for r agree 
with the exact result r = 5/2. Above the upper critical 
dimension, which is the case for our models, all critical 
exponents can be expressed in terms of a single one. So 
we arrive at the same critical exponents as in ordinary 
percolation. 


IV. SUSCEPTIBILITY AND ORDER 
PARAMETER 


According to Ref. 0 , the order parameter and gener¬ 
alized susceptibility for this class of problems are related 
to the probability c-i that two nodes selected by the model 
rules fall within the same cluster, 

C2 = E^f + [ 1 -( 1 - S n 2 -x/^V + 0 2 . (22) 


The first term on the right-hand side is the probability 
that both selected nodes belong to the same finite cluster, 
which is equal to the susceptibility \ divided by N. The 
second term is the probability that both nodes belong in 
the percolation cluster, which is equal to the square of the 
order parameter O. In the models under consideration, 
the susceptibility near the critical threshold t c is 

* = E ,S p[g) =rn 2 Y,sP(s)=rn 2 (s)p, (23) 

S ' ' S 

where we used Eq. ©, and (s)p is the first moment of 
the distribution P(s). For t > t c , summing both sides of 
Eq. © over s we get 

sc 

— = 2[1 - (1 - SH £ sQ(s) = 2 m 2 S(s)p. (24) 


Similarly, for t < t c , multiplying both sides of Eq. © by 
s and summing over s we obtain 

sQ(s ) 

Using Eqs. (|2Tjl . (|24]l . and (l25l) we find that the first 
moment of the distribution P(s) is symmetric below and 
above t c , namely (s)p = (2m 2 )~ 1 \t—t c \~ 1 . Then, the 
asymptotics of the susceptibility is independent of to, 

X=\\t-t c \-\ (26) 




The critical singularity of the order parameter 0 = 1 — 
(1 — S) m = mS is 

O = 87rm 3 * ao(t — t c ), (27) 

where we have used Eq. (ELD- Notice that, in contrast 
to Xi the critical amplitude of the order parameter O 
depends on ao and to. 


V. ESTIMATES 


Let us estimate P(s,f) by substituting the approxi¬ 
mated relation ([3]) into the evolution equation ([1|), 

^ ~ w 2 s E t)P(s— u i t) — 2 msP(s, t). (28) 
oz 

u= 1 

which is valid for large s and t close to t c . Let us rewrite 
the last equation in terms of the rescaled distribution 
P(s, t) = mP{s 1 1) and time t = mt 7 


dP(s,t ) 
dt 


s-l 

= s P(u, t)P(s—u, t ) — 2sP(s, t). 

U= 1 


(29) 


We assume that Eq. (BSD. being asymptotically exact 
near the critical point, describes approximately P(s,t) 
in the full range of cluster sizes and time. This equation 
coincides with the exact Eq. © for ordinary percolation 
(to = 1), which, for the initial condition P(s,0) = d Sj i, 
has the solution P(s,t c ) = r)s -3 / 2 at the criti¬ 

cal point t c = 1/2. This readily leads to the following 
estimates for t c and ao of our problem: 


t 


e 

c 


i 

2 TO ’ 


(30) 


and 


a 


e _ 

0 — 



(31) 


We also estimate the critical amplitude B of the perco¬ 
lation cluster relative size, S = B(t — t c ). Inserting ag 
into Eq. (1211) gives B e = 4, independently of to. Table |T] 
shows the numerical results for t c , a o, and B for different 
to, and compares them with the estimates of this section. 
Notice that our simple estimate produces surprisingly ac¬ 
curate results for to > 1. The estimate t e c is especially 
good, with an error of only 5 to 10%, while the estimate 
clq has a relative error about 2 times larger. 


VI. CONCLUSIONS 

In the present paper we have demonstrated that two 
types of the local optimization rule result in contrast¬ 
ing effects. The original Achlioptas rule based on selec¬ 
tion of the smallest clusters for merging together dras¬ 
tically changes the critical features of continuous phase 
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(a) (b) 




FIG. 2. Results of numerical solution of evolution equation © for s < 10 5 and different m. (a) Evolution of relative size of the 
percolation cluster S(t). For each m, the percolation cluster emerges with an exponent /3 = 1 at the critical point f c (m) given 
in Tabled] (b) Rescaled distribution P(s,t) shown in terms of aom 2 s 1/,2 (f — t c ). All curves P(s,t,m) for 14 combinations of 
s — 10 4 , 10 5 , and m = 1, 2, 3, 4, 5, 10, 20, collapse into one, exp[—47raom 4 s(t — t c ) 2 ], see Eq. O0l) . 


m 

T 

t c 

t e c = 1/(2m) 

do 

og = l/(y/2Trm) 

B 

i 

2.5000000(2) 

0.500000000(1) 

0.5 

0.398942(3) 

0.3989422... 

4.0000(1) 

2 

2.49998(5) 

0.2624198(1) 

0.25 

0.1755(1) 

0.1994711... 

3.096(4) 

3 

2.49999(5) 

0.1775290(1) 

0.1666... 

0.1128(1) 

0.1329807... 

2.878(5) 

4 

2.50000(5) 

0.1340937(1) 

0.125 

0.0830(1) 

0.0997355... 

2.770(7) 

5 

2.49999(5) 

0.1077242(1) 

0.1 

0.0657(1) 

0.0797884... 

2.712(9) 

10 

2.49999(5) 

0.05430571(5) 

0.05 

0.0322(1) 

0.0398942... 

2.61(2) 

20 

2.50001(5) 

0.02726238(2) 

0.025 

0.01595(5) 

0.0199471... 

2.56(2) 


TABLE I. Values of critical exponent r, critical point t c , and critical amplitudes ao and B for different m in the case of the 
evolution starting from isolated nodes. These values are obtained from the numerical solution of Eq. © for s < 10 5 by our 
method Q. For comparison, the table shows the estimates fg = l/(2m) and ag = l/(\/27r m) found in Sec. m The estimate 
for B is independent of m, B e = 4. 


transition compared to ordinary percolation and delays 
the transition. Inverting this rule and selecting the 
largest clusters for merging, we arrive at qualitatively the 
same critical behavior as for ordinary percolation though 
with a percolation threshold at much earlier times, see 
Fig. ©b). We have obtained the scaling functions and 
critical amplitudes for different to. Interestingly, the crit¬ 
ical point and critical amplitudes obtained numerically 
are very close to our simple analytical estimates taking 
into account only clusters of large sizes. We have in¬ 
dicated the order parameter and susceptibility in these 
problems and verified the Curie-Weiss law for the sus¬ 
ceptibility. 

In summar y, w e have applied the approach developed 
in Refs. 0, HlDl ilz| to processes generated by inverting 
Achlioptas rule and quantitatively described the perco¬ 
lation transition in these models. One could also mix 
the two rules, original and inverse, in the same model. 


For instance, at each step apply one of the two rules at 
random. Our results suggest that this is similar to the 
combination of the Achlioptas rule and the interconnec¬ 
tion of random nodes, which leads to the usual explosive 
percolation effects |30l - [33| . We based our conclusions on 
a set of models convenient for analytical treatment. We 
expect however that these conclusions are qualitatively 
valid for a much wider class of processes with inverse 
Achlioptas rules. 
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